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Abstract 

A model for nucleation of second phase at or around dislocation in a crystalline solid is considered. 
The model employs the Ginzburg-Landau theory of phase transition comprising the sextic term in 
order parameter {rf) in the Landau free energy. The ground state solution of the linearized time- 
independent Ginzburg-Landau equation has been derived, through which the spatial variation 
of the order parameter has been delineated. Moreover, a generic phase diagram indicating a 
tricritical behavior near and away from the dislocation is depicted. The relation between the 
classical nucleation theory and the Ginzburg-Landau approach has been discussed, for which the 
critical formation energy of nucleus is related to the maximal of the Landau potential energy. A 
numerical example illustrating the application of the model to the case of nucleation of hydrides 
in zirconium alloys is provided. 
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I. INTRODUCTION 



Nucleation of second phase in the vicinity of elastic defects such as dislocations occurs 
in many alloys 

HQ. 

For example, in Al-Zn-Mg alloys, dislocations not only induce and 
enhance nucleation and growth of the coherent Laves phase MgZn2 precipitates, but also 
produce a spatial precipitate size gradient around them 3|-|5| . Another example is formation 
of a new phase in ammonium bromide (NH^Br), namely (3 — > 7 phase transition, which is 
observed to occur in the vicinity of crystal dislocations [6|. In titanium and zirconium 
alloys, used in aerospace and nuclear industries, the presence of hydrogen leads to hydride 
formation (TiH x ., ZrH x .) close to and on dislocations, causing embrittlement of the alloy, 
thereby reducing its performance during service jz, 8]. 

Cahn {2I provided the first quantitative model for nucleation of second phase on disloca- 
tions in solids using classical nucleation theory. Calm's model assumes that a cross-section 
of the nucleus is circular, which is valid for a screw dislocation. Moreover, it posits that the 
nucleus is incoherent with the matrix. The issue of the formation of coherent nucleus on 
or near an edge dislocation has been studied theoretically by Lyubov & Solovyev 
Dollins 11] . These theoretical approaches have been thoroughly appraised by Larche 



fl and 



12|. In 



a recent study, Hin and coworkers [3] studied heterogeneous precipitation of FeC particles 
on dislocations in the iron-carbon binary system using kinetic Monte Carlo technique. 

Here, we present a generic model for nucleation of a new phase near edge dislocations in 
crystals. The model rests on the Ginzburg-Landau theory of phase transition in which an 
order parameter designates the symmetry of the system. Moreover, the elastic property of 
the solid is taken into consideration by the striction term in the free energy, which accounts 

nn 

for the interaction between the order parameter and deformation [14L [151. The model is in 

n n 

line with earlier approaches by Nabutovskii & Shapiro [16] and Boulbitch & Toledano |l]; 
where herein note the case of phase transition near edge dislocations has been elaborated. 
The model is pertinent to systems in which second phase nucleation is accompanied by a 
preferred orientation of nuclei under external force. This includes o/'-phase precipitation in 
Fe-N alloys {itI ]. #'-phase nucleation in Al-Cu alloys 18 1, 5-hydride formation in Zr-alloys 
containing hydrogen 19[. We should, however, point out that in this note, we only treat 
the details of the ordering (orientation) aspect of the problem, which is characterized by a 
non-conserved order parameter. That is, the effect of composition field is decoupled from 



the Ginzburg-Landau model. A more general formulation with coupled non-conserved order 



parameter, conserved variable (concentration) and elastic field was presented elsewhere 20], 
see also 



The paper is organized as follows. The model set up and the basic equations are described 
in section [TTl The ground state solution of the linearized steady-state Ginzburg-Landau 
equation, in the vicinity of edge dislocation, is derived in section IIHI using the method of 



Dubrovskii 2l|. The phase diagram ensued from the model is presented in section IIV1 
Section |V] discusses the relation between the present model and the classical nucleation 
theory. That section also includes a numerical example pertinent to nucleation of hydrides 
in zirconium alloys. Finally, in section IVII[ we end the paper with some concluding remarks. 
Some mathematical details are presented in the appendices. 

II. MODEL DESCRIPTION 

We consider an edge dislocation, the line of which coincides with the Oz axis and the 



Burgers vector with components b y = —b, b x = b z = 0, see 
crystal, such a dislocation creates a deformation potential 



igUJ In an elastically isotropic 



22] 



nr,0) = £— , (1) 
r 

where (r, 9) are the polar coordinates at the point of observation in a plane perpendicular to 
the dislocation line. Here B is a material dependent parameter denoting the strength of the 
potential. It is related to the basic constants of the metal, namely, B = bep(l — 2^)/[37r(l — 
z/)], where Ef is the Fermi energy and v is Poisson's ratio. 

The considered phase transformation is the nucleation of second phase in a solid solution. 
It is characterized by an order parameter accounting for the symmetry of the structure 
(system) under consideration. In general, it is defined by a vector field rj(r,t) being a 
function of space r and time t. We write the total free energy for the system 

T — J^st + F e l + J 7 inti (2) 

where T s t is the structural free energy, To, the elastic strain energy, and Tint is the interaction 
energy between the structural order parameter and the strain field. The structural free 
energy is 

= / [§(V»7) a + Vfo)]dr. (3) 



r sinO 




where the space integral is within the volume of the system. Here g(Vr]) 2 accounts for the 
spatial dependence of the order parameter, g is a positive constant, and V(rj) = ^r r] 2 + 
\uqTi a + ^v r] 6 is the Landau type potential energy, where r and u are, in general, functions 
of temperature and external field, here stress. The coefficient vq is taken to be a positive 
constant. Its role is to ensure stability. Also, we wrote rf = rj ■ rj and so on. The elastic free 
energy is 

T A = I [f (V ■ u) 2 + " ■ u) 2 ]dr, (4) 

where K and M are the bulk and shear modulus, respectively, U\j = (VjUi + ViUj)/2 is the 
strain tensor with V* = d/dxi, d the space dimensionality, and i,j stand for x, y, z in d = 3 
(x,y in d = 2). Finally, the interaction energy is 

J^int = K / rf V • u dr, (5) 



where rf V ■ u accounts for the interaction between deformation and the order parameter and 
the constant k, referred to as the striction factor, denotes the strength of this interaction. 
For structural orientation of the second phase, the order parameter may supposed to be a 
two-component vector field r\ = (j]x, r] 2 ), where (rji, r] 2 ) = (±//o> 0) would denote one preferred 
orientation of nuclei (precipitates) and (771,772) = (0, ±770) another, where 770 is a non-zero 



constant 



]. Furthermore, rj = (0, 0) would describe the solid solution with no precipitates 
present. Here, for the sake of simplicity, we assume that the order parameter is a scalar field 
(an Ising model), taking values of 77 = (solid solution) or 77 7^ (nucleus) 



The evolution of the order parameter is described by the time-dependent Ginzburg- 
Landau (GL) type equation 

where L a is the mobility; also the effect of the background random thermal (Langevin) noise 
is ignored. This is the basic kinetic equation for a non-conserved field, and it corresponds 



to the model A in the classification of Hohenberg & Halperin 23]. Moreover, defects in a 
crystal such as dislocations give rise to internal strains in the solid affecting the mechanical 
equilibrium condition. The equilibrium condition in the presence of an edge dislocation and 



the force field generated by the order parameter is expressed as (see Appendix |A~J and 241 ] ) 



MV 2 u + (A - M)VV • u + nS7rf = -Mbe y 5(x)5(y), (7) 

where A = K + 2M(1 — 1/d), b is the magnitude of the Burgers vector, e y denotes the 
unit vector along the y axis, and 8(») is the Dirac delta. Equation is used to solve u 
in terms of rf and the displacements of the dislocation; thereafter, we eliminate the elastic 
field u from the expression for the Euler-Lagrange condition on the free energy (Appendix 



lAjl . Hence, the total energy can be expressed as [20] 

/111" 
|(Vr/) 2 + -r lV 2 + -urf + -v oV 6 \dr, (8) 

where To = ^[0] is a function of temperature and stress, r\ = r — B cos9/\r\, r = 
a(T/T c — 1) = ar, a is a positive constant, T the temperature, T c the phase transition 
temperature in the absence of elastic coupling, u± = Uq — 2k, 2 /A for an edge dislocation 
embedded in the matrix, uq can be both positive and negative, and dr = rdrdOdz. Also, 
B = kA and A = (2b/7i)M/A = (b/ir)(l — 2z/)/(l — u). We regard J 7 [77] to have dimensions 
of energy [ML 2 T -2 ]. By considering 77 to have dimensions of inverse area [L~ 2 ], the variable 
ri gets dimensions [ML 3 T~ 2 ], and so on. Some of the constants appearing in the coefficients 
of the powers of 77, such as T c , the elasticity constants, the Burgers vector, are directly 
measurable for a particular system. Other parameters, such as g, a^, Vq, k, if not directly 
measurable, may be determined from ab initio type methods; see also section IVI1 
Let us first express Eq. flH]) in dimensionless form 

m = Fo J [{V^) 2 + \(U{p)-E)4, 2 + \^+ l -^]dp, (9) 



with dp = pdpd8d( and the introduced dimensionless parameters 

p = r/£, iP = - Vv ^/2B, C = z/£, U{p) = -cos 6/ p, (10) 

E = -r g/2B 2 , £ = g/2B, F = Bg/u h v = gul/8B 2 . (11) 

Here, we treat nucleation in the presence of the edge dislocation in thermal-mechanical 
equilibrium. Furthermore, we suppose that there is no Langevin noise in the system. Hence 
Eq. (JBJ), using Eq. reduces to 

V 2 ^= {U(p,8) -E)^ + 2^ + ^\ (12) 
2 d 2 Id Id 2 
^ dp 2 pdp~^ p 2 d0 2 ' ^ ^ 

and the boundary conditions: ip(p = oo, 6) = and ijj{p = 0, 9) ^ oo. 

One can show that for a certain E < 0, i.e. T > T c , Eq. ( 11 2 J) has a nontrivial solution 

for E > Eq (T < Tq). In the next section, we shall specify and calculate Tvo and obtain the 

bound state solutions of Eq. (I12p as — > Eq. 

III. BOUND STATES 

The value of E = Eq is found from Eq. ( f!2l) for which a solution would first emerge. This 
value is determined from the solution of the linearized Ginzburg-Landau equation, viz. 

V 2 ^=(U(p,9)-E)^ (14) 

which is the lowest "level" of Eq. ( ITS"]) , and it corresponds to the Gaussian approximation of 
the Ginzburg-Landau free energy functional. Equation f[T4"|) is equivalent to the Schrodinger 
equation in two dimensions. The bound state solution ip E to this equation exists if F^e] < 



and that occurs at E = Eq 25]. 



The variables in Eq. ( TT4|) with U(p, 6) defined in Eq. ( fTO]) do not separate. To solve Eq. 



( !T4|) we chose the method proposed by Dubrovskii 2l| outlined in Appendix |B] Recall that 
the eigenfunctions of the equation must satisfy the aforementioned boundary condition on 
ip. So we may write Eq. ( IB 8 1) in Appendix El expressed in terms of cosine-elliptic Mathieu 
function ce m (a, q), cf. Appendix [C] as 

Vwi = -^ 1 ce m (a,'i9,g)p At exp(-/3 ri p)iTi(-n,2 y u + l,2fi n p), (15) 

'71 



where A nm is the normalization constant, \x = \Ja + «i/2, d = {9 — 7r)/2, a is a separation 
constant, g, ai are variational constants to be determined, /3 n is a certain function of q, n is 
an integer, m is an even integer, and •, •) is the confluent hypergeometric function of 

the first kind. We should point out that for each q, the separation constant a = a m assumes 
an infinite set of discrete values depending on the parity and the index of the function 
ce m (a, q), cf. Appendix ICl 

The composite ground state eigenfunction corresponding to Eq. flT5|) is then 

ipoo(p,$,q) = ^ce (a i$,q)p tl exp(-/3p), (16) 

'IX 



where we put i-Fi(0, 2fi + 1, 2/3 p) = 1 and (3 = (3 Q . We should note that (3 and /i are functions 
of q, «i and «2 (see below). The normalization constant is evaluated to be 



,2(20)™* 

Ao °-^T(2 + 2^ (17) 

The ground state eigenvalue of the total Hamiltonian operator H = V 2 — U(p, to be 
minimized, is expressed as (E 00 <^ Eq) 

poo f* n /2 

E m ^ (00|#|00> = / pdp i; 00 (p,$)Hij 00 (p,$)d$. (18) 

J0 J-w/2 

Next, the integration of the right hand side of Eq. ( 1T81) over the variable $ yields 

(2B) 2+2fl f°° 

Eoo ^ f^T2p) J /+1 e M-Pp)W exp(-Bp)dp, (19) 

with t = _^ \d K^ Kl-a^q) 

dp 2 pdp Ap 2 2qp K J 



where ao(q) < is the eigenvalue of Eq. ( 1B3I) in Appendix (Bj corresponding to the eigen- 
function ceo(a, <?) and 

2 r' 2 , „ , d 2 

I 

-it/2 

Now, /C0(p) = E 00 (f)(p) with (p(p) = p M exp(— f3p) leads to 



2 f w/2 d 2 
K 2 = / ce (ao,^,q)—ce {a ,'d,q)d'd. (21) 

7T J-tt/2 



p= K o-^) fl= ^ i?oo = -/3 2 . (22) 
2g(AT + 1) 2 ' 00 1 K ' 



Additional analysis 



2l| yields (cf. Appendix] 



«i = K 2 - a , a 2 = ^i. (23) 



where «2 is an additional variational constant related to oq. Figure |2] depicts E Q0 (q) as 
a function of q. The obtained numerical values for the parameters, upon minimization of 
E 00 (q), are listed in Tabled! And in Figure EJ we depict the radial dependence of ^ oG°)$) 
at q = <?min- The eigenfunction ipoo{p, "$), corresponding to the normalized order parameter 
in the vicinity of the dislocation, is shown in Figure H] on the x?/-plane. 




FIG. 2. The ground state energy of the total Hamiltonian, Eq. (|19p . versus q. The minimum 
occurs at Q m i n = 3.918 and £ o(<?min) = -0.105, cf. Table [Q 



IV. PHASE EQUILIBRIA 

Let us study the phase diagram for the system under consideration. To first approxima- 
tion, as in [lj], the equilibrium field parameter ip is taken to be the ground state solution 
of linearized Eq. (|T2|) . i.e., ip = ^oo(p, Qm)- We then express the order parameter for the 

TABLE I. Numerical values for the parameters obtained by minimization of Eoo(q), Eq. (|18p . 
where q m = q^ 



limn- 



q m = 3.918 E 00 (q m ) = -0.105 K (q m ) = 1.284 

A 0Q (q m ) = 0.428 p(q m ) = 0.3245 a (q m ) = -4.159 

H(q m ) = 0.642 ai {q m ) = 5.808 a 2 (q m ) = 0.741 



S 



z=0 



z=n/A 



0.014^ 
0.012 
0.010 
0.008 
0.006 
0.004 
0.002 




z=n/2 





<Aoo 



0.014 
0.012 
0.010 
0.008 
0.006 
0.004 
0.002 




FIG. 3. The ground state eigenfunction ipoo(p, z, <? m in) with q m \ n = 3.918, cf. Eq. (fT6|) and Tabled 




FIG. 4. The distribution of the ground state eigenfunction ip = V'oo(P) Qmin) with (/ m i n = 3.918, 
cf. Eq. (fl~6|) and Table [U corresponding to the spatial variation of order parameter near the edge 
dislocation. Here x = pcosi? and y = psini?. 
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system in the form 



4' 

Sir 



rj(p, 0) = m^ce (a , q) ( — ) ' exp(-p/p ) + 0(m 



(24) 

'tt - \ Po ; 

where we used Eq. (fTB"]) . m is an amplitude, po = and A' Q0 = A 00 Pq. Note that if m is 
considered to be dimensionless, then the right-hand side of Eq. (124)) needs to scaled by the 
factor y/u\g/2B to make r] to have dimensions of inverse area [L -2 ]. 

Now substituting Eq. ( 1241) into (IE]), using the data in Table HJ and integrating over the 
volume (0<p<oo,0<$<7r, < z < L) yield the equilibrium free energy of an ordered 
nucleus around the dislocation. So the free energy is expressed in the form 

1 



jF(m) = T + 1 L 



4 b 



(25) 



where J-"o is the free energy of the defect in the parent phase (m = 0), L is the size of the 
crystal in the ^-direction, Uq = 2n 2 /A, and 



2B 



■C3 



OO P7T 



£ 2 ^^ ^ 2p7r ; J/i(l + 2p)' 
00 ceo(a ,'(9,g) — exp(-p/p ) 



2/i 



pdpdfi, 



(26) 
(27) 



'7T Vpo- 

where c 2 = 2.5903 and C3 = 1.1642 are obtained by appropriate integrations carried over the 
angle (0 < •& < it). Also the integrations in Eq. (|27|) can readily be evaluated, viz. 



1. 

0.1842 
0.0791 



4 1 " 2 ^ r(2 + 4p) 

pI r(2 + 2p) 2 ' 

2 7 x 3 _2-6m r ( 2 + 6p) 

— r(2 + 2p) 3 ' 



Pi 



(28) 
(29) 

(30) 



Equation ( 1251) is a kind of a mean field variant of the Landau free energy 26| with the 
order parameter m. The parameter v in Eq. fT25|) is taken to be a positive fixed constant, 
but f = (r — Tq) and u = (m — Mq) are varying parameters, assuming both positive and 
negative values. The system has a phase transition at f = 0, which can be either of first 
order or second order depending on sign of u. If u > 0, then tq > Tq gives m = 0, i.e. a 
situation with no nucleus. At tq = r^, nucleation occurs by second order transition; and 
for < ro < Tq, m 7^ 0, implying that the nuclei grow continuously. On the other hand, if 
u < 0, the transition is first order and the nuclei form on the coexistence line of the phase 
diagram, see below. 
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The value of m can be determined by the equilibrium condition ^(m) = J~o, and at the 
same time dF /dm = 0, which give 



rf = 4ffl(!!l4|. (31) 
4 (^ 6 ) V v 



Substituting this in J-"(m) — To, we find 



with Tq 



ro = r + - w -(n -n ), (32) 
--(l + ^)l^oo. (33) 



2/i 



/i V 27T/T 



This is the equation for first order phase transition line with roo denoting the ground state 
value of r . Next minimizing ^(m), i.e. dT /dm = and d 2 T /dm 2 > 0, then for r < r$ and 
Uq < Uq, we can calculate the resulting m (with admissible solutions); thereby leading to 

1 (ib 4 ) 2 

r ° = r "° + ^W"" ~ <Y ' (34) 

where the equality marks the line of instability, i.e. the onset of nucleation. Note that for 
ro < Tq and uq < Uq, one solution (m + ) gives two minima while another (m_) gives the 
maximum of the free energy. Also, (r — Tq) = const, and Uq > u* (and vq > 0) mark out 
the second order phase transition line. To illustrate this, we rewrite Eqs. ( 1321) and ( l34|) in 
a more concise form 

-u , rS , : - u . (35) 



Recalling that in the crystal bulk r = a(T/T c — 1), Eq. (!35l together with the line f = 
(or f = const.) and u > 0, define the phase diagram for the system in the (u, f) coordinates. 
Figure [5] shows such a diagram. Note that for a defect free crystal f = ro, while for a rigid 
crystal u = Uq. The point at which the first order transition turns to second order is the 
tricritical point designated by open circle in Fig. |5j Figure |6] shows a generic sextic Landau 
potential with coefficients u < and v > describing a first order phase transition between 
the parent-phase and the second-phase, denoted by I and II, respectively. 

V. CONNECTION TO CLASSICAL NUCLEATION THEORY 

Classical nucleation theory (CNT) has been used in the past to study the problem of the 



formation of coherent embryo on or near an edge dislocation 

11 



10 



12j | . In this framework, the 




FIG. 5. Phase diagram of nucleation in the bulk of rigid crystal (uQ,ro) and near the edge disloca- 
tion (u, f) obtained from Landau's potential (see figure[6j). The dashed line describes the first order 
transition, the solid line the border of instability, and the open circles are the tricritical points. 



free energy of the formation can be expressed as 12 1 



2tt Att 

AF = KAR 2 + 4vr 7j R 2 /i effJ R 3 , (36) 

3 3 

where R is the radius of the nucleus, K is the bulk modulus, 7 is the surface tension of the 
nucleus in the absence of dislocation, and /i c g- is the free energy difference per unit volume 
between the metastable and stable phases. The parameter A was defined earlier, i.e. after 
Eq. OH]). Now, AF takes a maximum at R = R c given by 

4ir / KA\ 9 „ . 

, ^ 2 / KA\ 
with R c = ( 7 -__). 38 

Note that A has dimension of length and 7 that of energy per unit area. It is seen that the 
dislocation simply shifts the surface tension to an effective surface tension 7 e g = (7 — KA/6); 
and so we may write AF C = (47r/3)7 e ffi? 2 and R c = 27 cff //x ef f. 

12 



FIG. 6. The Landau potential energy V = |ry 2 + |r/ 4 + |r^ 6 , with v > 0, u < 0, r > 0, and r\ 
scalar. The curves determine the lines plotted in figure [5j The dashed curve marks the emergence 
of metastable second-phase II (two local minima) in the parent phase I (global minimum); the 
dotted line indicates that phases I and II are equally stable; and the solid line refers to when phase 
I is less stable than phase II, which remains in that state until r < 0. 

The surface tension of the nucleus 7 is the excess energy stored in the matrix/nucleus 
interface region per unit area. It can be expressed in terms of the coefficients of the Ginzburg- 
Landau free energy, cf. Eq. ([8]). Let us for the sake of simplicity ignore the positive 
stabilizing term rj 6 and put A = 0. Moreover, consider tq < and a planar interface whose 
normal direction is in the x-direction. The Ginzburg-Landau equation becomes 

d 2 rj , 
g-r-z + \ro\r] - U\T) =0. (39) 
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Then we have a well-known solution for the interface profile 

^ )= © 1/2tanh ©' 



(40) 



where £ = (g/2|ro|) 1//2 is a characteristic length (the correlation length) for the formation of 
the new phase. The surface tension is given by 



7 



dx. 



Here Vrj = drj/dx and 



\ «1 ( 2 r 



(41) 



(42) 



where the free energy density — r^/Aui at x — ±oo has been deducted from the total energy 
density. Utilizing now Eq. (1T2"]) . after some manipulation, we can write 



/oo 
(Vrjfdx, 
-oo 



(43) 



Next, substituting the interface profile solution Eq. ( 1401 in ( )43|) . then evaluating the integral, 
we obtain 

9 IrvJ . l9r-J3/2 

(44) 



2 |r |. 1 1/2 |2r | 3 / 2 



3 «i 3ni 
Thus, we can express the maximal free energy of formation 



with R r 



47T 

AF C = — 
3 

2 

/^cff 



(45) 
(46) 



On substitution of Eq. ( 146]) in (145]) using the formula for £, the maximal energy of 
formation can be expressed in the form 



with q 



2 12 a/2tt GtVq) 1 / 2 

3 4 Mi 

lui 

4r ^Irof 



l-^) 3 ( 



3 / y max 



(47) 
(48) 



and V max = rQ/4w! is the maximum value for the Landau potential ( 1421) . Therefore, the 
maximal value of the formation energy, in the CNT setting, is proportional to the square of 
the maximum of V(rj). Note also that in the absence of defect g = 0. 
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VI. DISCUSSION 



Let us first discuss the characteristic length of the system under consideration. We recall 
the radial part of the bound state solution, Eq. (|16p . namely, 

(j)(p) = p M e -P/Po ; (49) 

where p = or more explicitly 

Po = /=== = ( ) ■ (50) 

To obtain Eq. (15 Op . we utilized the relations in ( I2"2"j) and ( ITTj) . with the ground state value 
of r , i.e. r o oc -Eqq 1 . Scaling p with £ to attain dimension of length, we get 

^po=(^v /2 . (si) 



■2|r 00 | 

The characteristic length £ ma Y be interpreted as the ground state size of the embryo; and 
since r = ar, we may write r 00 = a r , where r = T /T c — 1; thereby £ — (5 r /2« ) 1/ ' 2 |tq |~ 1//2 , 
with To being the ground-state nucleation temperature. Thus the size of the nucleus at the 
onset of nucleation is finite [l|. 

A point worth noting is that the Landau mean field theory of phase transition is valid 
so long as the fluctuations of the order parameter in a volume with linear dimension of 
order £ is small compared with the characteristic equilibrium value fj = {\tq\/ux) 1 ' 2 • The 
applicability of the mean field thermodynamics to describe the fluctuations at the onset of 
nucleation may be checked through the Levanyuk-Ginzburg criterion (see e.g., j^j]), which 
states 



« 1 7b | « 1. (52) 



The ratio on the left-hand side of Eq. (152]) is dimensionless (/cb is the Boltzmann constant) 
and is referred to as the Ginzburg number Gi. This number can be expressed in terms of 
characteristic lengths in the system, i.e. Gi = (£g/£o) 6 ) where 



So for the mean field theory to be valid £g << Co- We may also express the surface tension 
in terms of the Ginzburg number in the manner 

15 



Thus, in this formulation, 7 is meaningful so long as | to | >> Gi. 

Now to make the concepts more tangible, we give a numerical example relevant to pre- 



cipitation of hydrides in zirconium alloys with a terminal solid solubi 
of about 2 weight parts per million hydrogen at room temperature 



ity for precipitation 
. Typical values for 



Poisson's ratio v = 0.37, Young's modulus Y = 100 GPa [29|, and the magnitude of the 
Burgers vector b = 0.25 nm 30| can be used. Furthermore, taking k/Y = 0.1 m 4 , g = 10~ 8 
Jm 3 , ao/T c = 2.2 x 10 5 JmK -1 , we find To — T c « 10 K. With these numerical values, the 
critical radius for the formation of embryo nucleus, according to Eq. (IBTj) . turns out as 
£0 ~ 50 nm. Next, the surface tension of the nucleus can be estimated by using Eq. 
If Mi ~ 10 6 Jm 5 , then 7 rs 0.3 Jm -2 . Finally, with the aforementioned numerical values 
Gi = 2.6 x 10" 13 . 

For an application of the model to a particular material, such as precipitation of second 
phase ZrHa, in Zr matrix, the chemical free energy functional, expressed in terms of the 
concentration of hydrogen in the matrix, should be coupled to the system of the governing 
equations. In this context, an additional basic equation, namely the diffusion equation or the 
Cahn-Hilliard equation, should be solved together with the Ginzburg-Landau equation 20 1 . 
Moreover, the effect of appreciable plastic work done by the hydride under stress through 
volume expansion (Tr(e^) « 0.16) needs to be taken into account. Consequently, even more 
model parameters will appear in the equations, where their values need to be determined. 
As noted in section [Til some of the parameters, such as the elasticity constants, the Burgers 
vector, etc, are directly measurable. Other parameters, such as the coefficients of the Landau 
expansion may be calculated by ab initio type methods. For example, for the precipitation 
of 6' phase in defect free Al-Cu alloys such analysis has successfully been carried out [31]. 
For the aforementioned numerical example, the determination of the parameters ao, g, u±, 
and K would be worthwhile. 



VII. CONCLUDING REMARKS 

In this paper we have presented a model for nucleation of second phase at or around 
dislocation in a crystalline solid. The model employs the Ginzburg-Landau approach for 
phase transition comprising the sextic term in order parameter (77 s ) in the Landau potential 
energy. The asymptotically exact solution of linearized time-independent Ginzburg-Landau 

16 



equation has been found, through which the spatial variation of the order parameter in the 
vicinity of dislocation has been delineated. Moreover, a generic phase diagram indicating a 
tricritical behavior near and away from the dislocation is depicted. The relation between the 
classical nucleation theory and the Ginzburg-Landau approach has been evaluated, for which 
the critical formation energy of nucleus is related to the maximal of the Landau potential 
energy. In a numerical example, a number of model parameters is fixed to certain numerical 
values to obtain plausible results, namely, the shift in the phase transition temperature due 
to the presence of the dislocation, the critical radius for the formation of embryo nucleus, 
and the surface tension of the nucleus. The numerical values for the model input constants 
need to be determined by detailed experiments and/or ab initio computations for the system 
under consideration. 

Numerical calculations applicable to the time-dependent Ginzburg-Landau equation with 
the nonlinear terms will be presented elsewhere. Subsequent steps in our study are an 
extension of the calculations to a two-component field structural order parameter, the space- 
time evolution of the order parameter in the vicinity of defects, and the coupling of the 
composition and the structure order parameters. 
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Appendix A: Elimination of elastic field 

The procedure of eliminating the elastic field appearing in the expressions of the system 
free energy has been discussed by many authors in the literature, e.g. (3, 32, 33 1. Here, we 



outline a simple procedure applicable to our case. The total free energy functional, Eq. (j2J), 
is expressed in terms of two field variables, the order parameter r\ and the displacement field 
vector u. In equilibrium, the spatial distribution of the order parameter and the displacement 
are determined by the Euler-Lagrange equation: 

£ = < A1 > 

b -f = Qi, (A2) 
17 



where Qi on denotes the force field due to the presence of an elastic defect in the solid. It 



has been shown in 



24( that this force per unit volume element dV is Q/dV = Mr x b<5(£), 



where £ is a two-dimensional radius vector perpendicular to r with origin at the dislocation 
line, r is a vector parallel to the dislocation line, b is the Burgers vector, and S(») is the 
Dirac delta. Thus, by virtue of Eqs. ©-([S]), Eqs. (lAip and (1A1I) . respectively, yield 



r ?7 + Uoi] + v r] + 2m] V • u = gV rj, 
MV 2 u + (A - M)VV • u + kVt] 2 = Mr x b<5(£), 



(A3) 
(A4) 



where A = K + 2M(1 — 1/d). For an edge dislocation, r = — e z is constant along the 
dislocation line, while the Burgers vector is in the x direction b = e x . Here, denotes the 
unit vector along the i axis. Hence the right hand side of Eq. (IA4h for an edge dislocation 
becomes — Mbe y 5(x)S(y), which gives Eq. (JTj). Solving now Eq. ([7j) the two components of 
the displacement vector are 



u, 



b 

2^ 
b 

2^ 



D 



x 2 + y 1 

y 2 



D - - - — arctan ( — 

M rx 2 + y 2 
x- + y 2 ~ ~ 2A ° g V 52 



2„2 



A 



2„2 



(A5) 
(A6) 



where D = (K + M/3)/A and V 2 is the inverse Laplacian operator. The dilatation strain 
V • u = V x u x + Vyiiy is (cf. with relations in 24j) 



Vu 



b 2M / y 
2tx A \x 2 + y 2 J A 



Inserting Eq. (IA7I) into Eq. flA3[) . we write 



?7 • 



(A7) 



gV 2 rj = r\T\ + uirf + v ?T 



(A8) 



where renormalized parameters are 



7*1 = Tq — Kv4 



cos 



Ml = M 



2^2 

"X' 



(A9) 



with A = (2b/ir)M/ A, 9 = arctan(|//x) and r 2 = x 2 + y 2 
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Appendix B: Eigenfunctions and eigenvalues of H: Variational and perturbative 
method 



We use the method proposed by Dubrovskii [21] to solve Eq. (|T0|) . In this method one 
first writes Eq. (fT4"j) in the operator form Hip = Eip, where H = H + Hi, and 

2 «i — 2qcos8 «2 



4p 2 p 
sin9 «! — 2qcos8 a 2 

P 



(Bl) 
(B2) 



p Ap 2 

where q, ot\ and a 2 are variational parameters to be determined. Now the variables in the 
equation H ip = E ip are separable. Hence we write ip(p,9) = i?(p)$(6 ) ), then the following 
equations are obtained 

d 2 $ , 



dz 2 

d 2 R 1 dR 
dp 2 p dp 



En 



- 2gcos2z)$ = 0, 

a 2 \ I 2 
P 



p 2 ± 



R = 0. 



(B3) 
(B4) 



Here, 9 = 2z + ir, I 2 = (a + «i)/4 and a is a separation constant. Equation (IB 3 1) is Mathieu's 
equation 34j,[35|]; whereas E q. (B4p is the 2-dimensional radial Schrodinger equation, and its 



solution depends only on |/| 36|. The periodicity condition of 9 with a period 2ir is satisfied 
by the Mathieu functions ce m (z,q) and se m (z,q) where m is the order of the functions 
assuming even integers here; hence solution to Eq. flB3j) can be expressed as 



0,2,4,... 



2,4,6, 



(B5) 



ce m (a,z,q), m 
se m (a,z,q), m 

where the notation ce and se comes from cosine-elliptic and sine-elliptic, respectively. One 
should note that for each q, the separation constant a assumes an infinite set of discrete 
values depending on the index and the parity of the functions; see Appendix O 

The normalized eigenfunction of the radial equation (1B4|) is well known, e.g. ref. 
and can be expressed in the form 

1 1/2 



Rnl(p) 



2f3 n 



(n+\l\ - 1)! 



(2n - l)(n - |/| - 1)! 



(2f3 n pfexp(-(3 n p) 



x J F 1 (-n+\l\ + 1,2111 + 1,2^), 



(B6) 



where 2(3 n = a 2 / {n — 1/2), n = 1, 2, 3, ... , and iF\(a, b, z) is the confluent hypergeometric 



function of the first kind 



341 ] . Also in this setting, / = 0, 1, 2, 3, 



n 



1. The bound state 
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energy levels are 



£l 0) 



4(ra- 1/2) 2 

Thus the composite eigenfunction for the system under consideration is 



mnl 



$ m (z,q)R nl (p). 



(B7) 



(B8) 



In general the eigenf unctions and the eigenvalues of the Hamiltonian Hq, see Eq. ( IBID , 
can be classified by three characteristic numbers {n,p, m}, where n = 0, 1, 2, . . . is the radial 
solution index, p — 0, 1 the parity index, determining the inversion symmetry with respect 
to variable z, and m is the index of the Mathieu function, taking up only even integers, viz., 
m = 0, 2, 4, . . . for p = and m = 2, 4, 6, . . . for p = 1. The eigenf unctions of H can be 
expressed in the form 



where $ r 



exp(-/3 npm p)iF 1 (-n, 2p + 1, 2p npm p)$ npi 



■ npm 



P 



npm 
/' 

a 

4 

and £(° p ) m 



'npm 

Op,0CG m -\~ 5p t \SC m i 

2/i + 2n + 1 ' 
^(a + « 1 ) 1/2 , 



V2(2^J 1+ ^ / r(2/x + 7i + l) 
r(2/x + l) V n! (2/i + n + l) 



"/3 



npmi 



(B9) 
(BIO) 
(Bll) 

(B12) 
(B13) 

(B14) 

(B15) 



where define the energy eigenvalues for the Hamiltonian (2l|, |25|]. Considering the 

ground state -^ooo with (3 = /3 oo = otij (2/x + 1) and /i = (ao + «i) 1 / 2 /2, and comparing with 
the relations in Eq. ( 1221) . the relations in Eq. ( 1231) follow. 

Next, using standard perturbation technique, the first (non-zero) shift is calculated 
through flBT5|) 

,W_ P (o)/^, ((000|^|002)) 2 



000 



with 



(000|F|002) 



-^ooo (lo 



pdp 



(0) 



tt/2 



V>ooo(<?o, p, z)Hip 002 (q2, P, z)dz. 



(B16) 



(B17) 



Here, g an d q 2 are the the minima of E^ (q) and E^ 2 (q), respectively, calculated through 
Eq. (IB15j) with numerical values listed below 
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State 






a 2 


J-J npm ^yminj 


|000) 


3.91792 


5.80765 


0.741165 


-0.105288 


|002) 


2.76301 


3.77263 


0.682702 


-0.0277153 



Hence, Eq. flETfj) gives (000|if|002) = 0.0221 and Eq. (1BT6]) yields -0.112. Our 

calculated numerical values agree well with those reported by Dubrovskii |21[ for the ground 
state, but differ somewhat for the first excited state. 



Appendix C: Angular Solutions 

The angular solutions § m (z,q), Eq. ( 1B5j) . expressed in terms of the Mathieu functions 
ce m (a,q, z) and se m (a,z,q) have subtle properties. The value of a in Eq. flB3j) form the 
eigenvalues of the equation; for ce m they are commonly denoted by a m (q), whereas for se m 
by b m (q). Hence we write: a =>- a(p, m) = 5 P)0 a m + 5 p ^b m with p — 0, 1 and 5 P)Q denoting the 
Kronecker 5. The index m takes up even integers starting from zero for p = and from 2 for 
p — 1. The Mathieu functions form a complete set of functions in the interval < z < 2tt, 
namely 

/•27T /*27T 

/ ce m ce n dz= / se m se n dz = n5 m , n (CI) 
Jo Jo 

Finally, since the Mathieu functions ce m and se m are periodic, they can be expanded by 
Fourier series 

oo 

ce 2r (a, q, z) = ^ -4. 2 fc(<?) cos(2/c2;) (C2) 

fc=0 

oo 

se 2r . +2 (a, q, z) = } -B 2fc+2 (g) sin[(2A; + 2)z\ (C3) 

fc=0 

where r = 0,1,2,.... The recursion relations between the coefficients are obtained by 
substituting these series into Eq. (1B3I) . For example, for ce 2r , we obtain 

qA 2 = aA (C4) 
q A 4 = {a- A)A 2 - 2qA (C5) 
g^2/t+2 = (a - 4:k 2 )A 2k - qA 2k -2 (C6) 

with k > 2. More on Mathieu functions is found in I34l. \3 
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